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The blackbody theory is revisited in the case of thermal electromagnetic helds inside uniaxial 
anisotropic media in thermal equilibrium with a heat bath. When these media are hyperbolic, we 
show that the spectral energy density of these fields radically differs from that predicted by Planck’s 
blackbody theory. We demonstrate that the maximum of their spectral energy density is shifted 
towards frequencies smaller than Wien’s frequency making these media apparently colder. Finally, 
we derive Stefan-Boltzmann’s law for hyperbolic media which becomes a quadratic function of the 
heat bath temperature. 

PACS numbers: 44.05.+e, 12.20.-m, 44.40.+a, 78.67.-n 


In 1901 Planck [T] derived the famous law describing 
the spectral distribution of energy of a blackbody (BB) 
by introducing the concept of quantum of light laying so 
the foundation of quantum physics. In his description 
of the problem, the electromagnetic field inside a cavity 
made with opaque walls which is set at a constant tem- 
perature is studied. In this formulation [5], the cavity is 
at thermal equilibrium and acts as a heat bath. The walls 
of the cavity emit and absorb electromagnetic waves so 
that the field itself becomes equilibrated. The internal 
energy density of the electromagnetic field in the cav- 
ity with volume V for both principal polarization states 
(abbreviated by s and p) is then given by 
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where h, and c denote Planck’s constant, Boltzmann’s 
constant and the velocity of light in vacuum, while P and 
C are Riemann’s gamma and zeta functions. Here and in 
the following we neglect vacuum fluctuations. 

In this Letter, we revisit this old problem when the 
cavity is filled with a uniaxial medium [3] with a relative 
permittivity tensor of the form 

/e+ 0 0\ 

e = 0 0 . (2) 

VO 0 ej 


modes (EMs) exist and satisfy the dispersion relations [5] 


' 9 ’ 

(OM) 

( 3 ) 

e_L e_L 



kl 

2 ’ 

^711 ^ 1 

(EM) 

( 4 ) 


where k± (A:||) is the wave number component perpen- 
dicular (parallel) to the optical axis. In usual dielec- 
tric uniaxial media the principal constants e± and e|| are 
both positive and the iso-frequency surfaces defined by 
relations ^ and ffl are spheres or ellipsoids, resp., as 
illustrated in Fig. [ifa). On the other hand, when e|| < 0 
and e_L > 0 or e|| > 0 and e_L < 0 the iso-frequency 
surfaces of the EM are two- or one-sheeted hyperboloids 
[see Fig. Hb)]. The first class of such uniaxial medium 
is called hyperbolic of type I while the second one hy- 
perbolic of type II Of course both e_L and e|| can 



Here without loss of generality we assume that the op- 
tical axis points into the z-direction; e|| is the permit- 
tivity along the optical axis and e_L is the permittivity 
perpendicular to the optical axis. For convenience, if 
not specified differently we neglect dispersion, dissipation 
and nonlocal effects in the following. Within such mate- 
rials so-called ordinary modes (OMs) and extraordinary 


Figure 1: (a) Cavity at temperature T containing an isotropic 
medium of permittivity e > 0 or an anisotropic (uniaxial) 
medium with e_L > 0 and e|| > 0. The particular case 
e± = e|| = 1 corresponds to the classical BB. (b) Cavity 
at temperature T containing a hyperbolic medium of type I 
(ej_ > 0 and e|| < 0) and of type II (e^ < 0 and e|| > 0). 
The isofrequency curves are plotted inside the cavities in the 
plane {k±, A:||). 
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also be negative. In such uniaxial metallic-like materials 
no propagating modes exist and the upcoming quantities 
are all 0. 

To start, we focus our attention on the electromag- 
netic field inside a cavity filled with an arbitrary uniaxial 
medium. The spectral density of states (DOS) defined 
as the energy density U normalized to the mean energy 
of a harmonic oscillator and associated to the thermal 
field can be calculated either by counting the modes in 
the wavevector space using expressions (3) and (4) or by 
means of the generalized trace formula mnim 

D{uj) = 4^ImTrreG™(r,r,a;)-f/xG™(r,r,w)l, (5) 

where and are the electric and magnetic 

Green’s dyadics for the bulk material and fi is the rela- 
tive permeability tensor. The result for a classical BB 
can be retrieved by using e = /x = 1 with the unit 
dyad 1. Then the above expression reduces to the well- 
known expression Dg’^^(w) = In the following, for 

the sake of clarity, we consider non-magnetic materials 
(i.e. /2 = 1) only. By inserting the general expres- 

sion of dyadic Green’s tensors of uniaxial materials [8] 
in the trace formula © it is straight forward to derive 
the DOS for the three different classes of uniaxial media. 
Assuming that those media are lossless then in dielectric 
anisotropic media the DOS for the OMs and Dfy for 
the EMs are given by the following expressions ISIEQ] 
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On the other hand, in the hyperbolic case we obtain 
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the EMs of type I and type II hyperbolic media coincides 
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With the help of the DOS we can determine the ther- 
modynamic potentials of the photon gas inside the uni- 
axial material. By definition, the internal and the free 
energy per unit volume are given by m 
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Einally, from the internal and free energy we can also 
determine the entropy per unit volume by 
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Clearly, by means of these expressions we can derive any 
thermodynamic property of the photon gas inside the 
cavity as the pressure P = ^, the photonic heat capacity 
Cv = fy, etc. 

Let us first have a look at the expressions for the ordi- 
nary uniaxial material. In this case we obtain 
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for the type I and type II media, respectively. Note that, 
we have introduced a cutoff wavenumber fc-L,max which 
for dispersive media can be a function of the frequency 
and which is determined by the concrete (atomic or meta) 
structure of the medium. Eor an ideal hyperbolic mate- 
rial fc-L,max is infinity so that the DOS diverges as was 
pointed out previously m- However, for any real struc- 
ture fc-L,max is a finite quantity m- For artificial hyper- 
bolic structures it is mainly determined by the unit-cell 
size of the meta structure. Note further that the DOS of 


Therefore, when e_L = e|| = 1 we recover the classical 
blackbordy result. And the relation between the internal 
energy, the free energy and the entropy have the familiar 
forms 
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Note that these relations are the same as for a usual BB 
because the DOS of the field inside a dielectric uniaxial 
medium is proportional to 

On the contrary, in type I and type II hyperbolic media 
we have seen that the DOS of the EMs is linear in ui as in 
a 2D photon gas in vacuum. It follows that the relations 
between the thermodynamic properties of the photon gas 
are radically different in that case. Indeed, we obtain 
(fc_L,max > y \/k|| I) 
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Hence U, F and S are proportional to and not any- 
more to T^. This result is a direct consequence of the 
linear dependence of the electromagnetic DOS inside hy- 
perbolic media with respect to ui. Naturally, for the OMs 
we find 


C/f = and C/fi = 0. (21) 


Note that for the type II hyperbolic material the inter- 
nal energy of the OMs is zero, since there are no OMs 
in such a material. The internal energy of the OMs in a 
type I hyperbolic materials is just the same as in a di- 
electric uniaxial medium. Hence, the relations between 
the thermodynamic potentials are the same as in a di- 
electric uniaxial medium. However, in typical hyperbolic 
(meta)materials the maximal wave vector is much larger 
than the vacuum wave vector fc_L,max “ making the 
material properties dominated by the EMs. 

Another consequence of the linearity of the DOS with 
respect to uj inside a hyperbolic medium is the spectral 
shift of Wien’s frequency Wmax (resp. wavelength Amax) 
at which the energy distribution function has its maxi- 
mum. For both type I and type H hyperbolic media we 
find after a straight forward calculation from relations 
(9) and (11) (fc_L,max “ Y^|e||I) that this maximum is 
reached when 
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whereas for a usual BB = 2.82 and 27r4/Amax = 

4.965. Here we have introduced the thermal coherence 
length lc = [5]. Hence we see that Wien’s frequency 

is shifted toward smaller values (i.e. the medium appears 
colder than a classical BB) and the maximum vacuum 
wavelength to larger values (see Fig. [^. 

It is now interesting to compare the internal energy of 
the EMs in a hyperbolic material with that of a classical 
BB. From expressions 0 and ([T^ we immediately get 
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If A denotes the unit-cell size of our hyperbolic material 
then fc_L,max = (27r)/A, so that 
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At a temperature of 300 K the coherence length is U = 
7.6 /im. The period of realistic artificial hyperbolic meta- 
materials is typically larger than A « 10 nm. In natural 
hyperbolic materials the unit cell size reduces to the in- 
teratomic spacing, i.e. A « 1 A. Hence, the internal en- 
ergy of thermal radiation inside a hyperbolic cavity can 
be 3 to 5 orders of magnitude larger than that of a perfect 
BB. The same is of course also true for the free energy 
and the entropy. This result suggests that the radiative 
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Figure 2: Spectral energy density (in log-log scale) of a type 
I and type II hyperbolic BB (red solid line) with |£||e±| = 1 
and fcn.max = ^ with A = 100 nm at T = 300 K. This dis- 
tribution is compared with the classical BB spectrum (blue 
dashed line). The solid and dashed straight lines show the 
asymptotic behavior in ca and of the hyperbolic and classi- 
cal BB spectrum. The arrows indicate Wien’s frequencies in 
both cases. 


heat flux inside a hyperbolic material is dramatically en- 
hanced compared to that of a classical BB. 

In order to evaluate the flux radiated by a cavity filled 
with a hyperbolic medium into a hyperbolic medium and 
to derive Stefan Boltzmann’s law we calculate now the 
Poynting vector in the cavity in the direction of the prin- 
cipal optical axis by assuming, for convenience, that the 
cavity opening (see Fig. [^ is along this axis. Using the 
framework of fluctuational electrodynamics theory the 
ensemble average of the Poynting vector (for any dis- 
persive and dissipative anisotropic medium) reads [H] 
(Einstein’s convention) 
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Here we have introduced the Levi-Civita tensor U/3'y 
the permeability of vacuum /iq. Note that (G^®(r,r') = 
X G™(r, r') and that ^ denote the conjugate trans- 
pose. In order to determine the heat flux, we assume that 
the cavity is infinitely large so that we can replace it by 
a uniaxial halfspace at a given temperature T. Inserting 
the Green’s dyadic [5] and integrating over this halfspace 
with volume V we find after a lengthy calculation (see 
Ref. [3] ) in the lossless limit, the relatively simple ex- 
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pression 
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for the mean Poynting vector or heat flux along the sur- 
face normal with 
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Evaluating this expression for the mean Poynting vec- 
tor for the dielectric uniaxial material, flrst, we have 

$i/ dwZ^(w,T)^--{^^}. (28) 

For non-dispersive materials this simplifles to 
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When e_L = e|| = 1 we flnd again the usual BB result, 
i.e. Stefan-Boltzmann’s law. On the other hand, inside 
a uniaxial material (as inside an isotropic material) with 
e_L > 1 and e|| >1 the radiative heat flux is larger than 
the BB value, which is a well-known fact [H] . 

In the case of hyperbolic media these results radically 
change. Before seeing this, let us first consider the OMs. 
For we find of course the same relation as for the 
dielectric anisotropic material, whereas as a consequence 
that there do not exist any OMs in a type II hyperbolic 
material we find $jj = 0. On the contrary, for the EMs 
we find 

nOO 1 

$f = du;U{u;,T)^^f^ (30) 

and 
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Hence, in the non-dispersive case, where fcn,max ^ 
yvWl have 

« $ej « ^|^^^^^(fcBT)2. (32) 

In this case, we see that the heat flux is proportional to 
and not anymore to T^ as in “classical” Stefan Boltz- 
mann’s law. Comparing this quantity with the classical 
BB results, we find 
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Hence the normalized heat flux is proportional to 
(fc_L,max^c)^ which is due to the fact that the heat flux 
scales like the area of the projection of the isofrequency 


surface in k-space [miis] or like the number of transver- 
sal modes [171 [18], resp. We have seen before in Eq. (23) 


that the ratio of the energy density of a hyperbolic ma- 
terial and that of a BB is only linear in fc_L,max4- This is 
quite astonishing, since for artificial hyperbolic materials 
with a unit cell size A of 10 nm and for natural hyperbolic 
materials with a unit-cell size A « 1 A we can now expect 
a hyperbolic BB heat flux 6 to 10 orders of magnitude 
larger than that of a usual BB at T = 300K. At cryogenic 
temperatures A becomes very large so that this ratio can 
become even much larger. To substantiate this state- 
ment let us consider a numerical example: We consider 
an artificial hyperbolic medium made of layered periodic 
structure of GaN and SiO^ slabs of thickness 5 nm each. 
Then, the period is A = 10 nm so that the theoretical up- 
per limit for the wavenumber is fc_L,max ~ 6.28 x 10®m“^. 
If we consider now two blackbodies made of and sepa- 
rated by such a material with Tmean = 300 K and small 
temperature difference AT then the transferred energy 
per unit area and temperature (heat transfer coefficient) 
for the EMs is A$®/jj/AT « 8.93 x lO^W/m^K. Note 
that this value for A<I>j^jj/AT has to be considered as 
an upper limit. It is interesting to compare this radiative 
heat flux to the heat conduction by phonons and electrons 
inside the hyperbolic multilayer structure. At ambiant 
temperature, the thermal conductivity of each unit layer, 
is KGaN < 0.5W/(mK) and KsiOa ~ 0.4W/(mK) POIET] 
the effective thermal conductivity of the whole struc- 
ture is about 0.44 W/(mK) when assuming that the ther- 
mal resistance of the combined multilayer structure is 
the averaged sum of the resistances of both materials 
(i.e. Kapitza resistances are neglected). Hence the heat 
transfer coefficient by heat conduction through a hyper- 
bolic multilayer structure of 200 nm (20 periods) is about 
2.2 X 10®W/m^K. Therefore, a hyperbolic BB can the- 
oretically have a radiative heat flux even larger than 
heat conduction in weakly conducting composite struc- 
tures [5^ HSj. It is important to note here that when 
a uniaxial BB radiates into vacuum the maxi- 

mum wave number is k± = w/c. Therefore for dielectric 
uniaxial materials we find, neglecting reflections at the 
cavity opening, again (if |e_L| > 1) the usual BB result 
= $ 33 , whereas for hyperbolic materials we find in 
this case 


4.//® = 4>3'']P and = 0. (34) 

The type I hyperbolic BB behaves like a perfect BB and 
the type H hyperbolic BB behaves like a perfect metal or 
a“white” body. 

To summarize, we have extended the BB theory to 
arbitrary uniaxial materials. For dielectric anisotropic 
media we have seen that the thermodynamic properties 
of the photon gas inside such media are very similar to 
that of a classical BB. On the other hand, when these 
media are hyperbolic, the spectral energy distibution 
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of radiation is shifted towards frequencies smaller than 
Wien’s frequency making these media apparently colder. 
We have also shown that in contrast to Stefan Boltz- 
mann’s law, the heat flux radiated by these media de- 
pends quadratically on their temperature. Nevertheless, 
the magnitude of heat flux carried by these media can be 
several orders of magnitude larger than the flux radiated 
by a classical BB and may even exceed the heat flux car- 
ried by conduction in superlattices. Detailed derivations 
of the above relations and the underlying assumptions 
as well as more detailed discussions will be given else- 
where 0. 
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